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Abstract. Detections of CO, CS, SO, C 2 H, HCO+, HCN, 
HNC, H2CO, and C3H2 are reported from LIRS 36, a star- 
forming region in the Small Magellanic Cloud. C 18 0, NO, 
CH3OH, and most notably CN have not been detected, 
while the rare isotopes 13 CO and, tentatively, C 34 S are 
seen. This is so far the most extensive molecular multi- 
line study of an interstellar medium with a heavy element 
depletion exceeding a factor of four. 

The X = A(H 2 )//co conversion factor is ~ 4.8 x 
10 21 cm -2 (Kkms -1 ) -1 , slightly larger than the local 
Galactic disk value. The CO (1-0) beam averaged column 
density then becomes Y(H 2 ) « 3.7 x 10 21 cm -2 and the 
density n(H 2 ) « 100 cm~ 3 . A comparison with X-values 
from Rubio et al. (1993a) shows that on small scales (R ~ 
10 pc) A-values are more similar to Galactic disk values 
than previously anticipated, favoring a neutral interstellar 
medium of predominantly molecular nature in the cores. 
The J( 13 CO)//(C 18 0) line intensity ratio indicates an un- 
derabundance of 12 C ls O relative to 13 C 16 w.r.t. Galactic 
clouds. 7(HCO+)//(HCN) and /(HCN) //(HNC) line in- 
tensity ratios are > 1 and trace a warm (Tki n > 10 K) 
molecular gas exposed to a high ionizing flux. Detections 
of the CS J=2-l, 3-2, and 5-4 lines imply the presence 
of a high density core with n(H 2 ) = 10 5 — 10 7 cm -3 . In 
contrast to star- forming regions in the LMC, the CN 1-0 
line is substantially weaker than the corresponding ground 
rotational transitions of HCN, HNC, and CS. CO, CS, 
HCO + , and H 2 CO fractional abundances are a factor ss 10 
smaller than corresponding values in Galactic disk clouds. 
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einmann@biaa21.biaa.sinica.edu.tw 

* Based on observations with the Swedish-ESO Submillime- 
ter Telescope (SEST) at the European Southern Observatory 
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Fractional abundances of HCN, HNC, and likely CN are 
even two orders of magnitude below their 'normal', Galac- 
tic disk values. The CN/CS abundance ratio is ^ 1. Based 
on chemical model calculations, we suggest that this is be- 
cause of the small metallicity of the SMC, which affects 
the destruction of CN but not CS, and because of the high 
molecular core density which also favors CN destruction. 

Key words: ISM: abundances - ISM: molecules - Galax- 
ies: abundances - Galaxies: ISM - Magellanic Clouds - 
Radio lines: ISM 



1. Introduction 

Recent studies of the dense molecular gas in five star- 
forming regions of the Large Magellanic Cloud (LMC) re- 
vealed a number of striking differences with respect to 
properties typically observed in the clouds of the Galactic 
disk (see Johansson et al. 1994; Chin et al. 1997): C 18 is 
underabundant relative to 13 CO, the C 18 0/C 17 ratio (« 
2) appears to be smaller than the canonical value of 3.5 in 
the Galactic interstellar medium, and HCO + to HCN line 
intensity ratios are larger than those in the Milky Way. 
In view of these results, one may expect that a thorough 
study of the Small Magellanic Cloud (SMC) with its low 
metallicity and strong UV radiation field will reveal even 
more drastic deviations from typical Galactic disk prop- 
erties, thus permitting insights into otherwise not directly 
discernible astrophysical and astrochcmical processes (cf. 
Johansson 1991, 1997; Rubio et al. 1993a, b, 1997; Lequeux 
et al. 1994; Chin et al. 1997; Heikkila et al. 1997). Em- 
ploying a Schottky receiver, Chin et al. (1997) reported 
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the detection of C2H, CS, and a tentative detection of 
CN toward LIRS 36, the source with strongest CO J— 1-0 
line temperature among investigated SMC IRAS sources 
(Rubio et al. 1993b). 

Here, a high sensitivity study of LIRS 36 is presented, 
providing a detailed molecular view of a star formation re- 
gion in the extremely metal poor environment (e. g. West- 
erlund 1990) that characterizes the SMC. 
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Fig. 1. Contour maps of LIRS36. a) 12 CO(1-0): contour levels 
are 0.0 to 9.0 by 1.5Kkms -1 integrated between dlsr = 121 
and 131 kms -1 . b) 13 CO(1-0): contour levels are 0.0 to 1.0 by 
0.2 K kms -1 integrated between vlsr = 121 and 131 kms -1 . 
c) 12 CO(2-l): contour levels are 2.5 to 17.5 by 2.5Kkms -1 
integrated between dlsr = 119 and 130 kms -1 . d) 12 CO(3-2): 
contour levels are 5 to 17 by 3Kkms -1 integrated between 
vlsr = 121 and 131 kms" 1 . Typical r.m.s. values are 0.14, 
0.25, and 0.3 K kms -1 , respectively. 



2. Observations 

The data were taken in January 1996 and 1997 and in July 
1997, using the 15-m Swedish-ESO Submillimetre Tele- 
scope (SEST) at La Silla, Chile. Two SIS receivers, one at 
A = 3 mm and one at 2 mm, yielded overall system tem- 
peratures, including sky noise, of order T sys = 250 K on 
a main beam brightness temperature (Tmb) scale while 
for 12 CO(1-0) T sys reached 400 K. In January 1997, a A 
= 1.3 mm SIS receiver yielded overall Tmb system tem- 
peratures of 500 - 1600 K at 218 - 245 GHz and 3000 
- 4000 K at 265 - 267 GHz in clear but humid weather. 



was employed. The backend was an acousto-optical spec- 
trometer (AOS) which was split into 2 x 1000 contiguous 
channels for simultaneous A ~ 3 and 2 mm observations. 
At A 1.3 and 0.85 mm, all 2000 channels were used to 
cover a similar velocity range. The channel separation of 
43 kHz corresponds to 0.04, 0.05 - 0.06, 0.08 - 0.10 and 
0.11 - 0.15 kms -1 for the frequency intervals 340 - 350, 
265 - 220, 150 - 130 and 115 - 85 GHz, respectively. De- 
pending on line frequency, the antenna beamwidth varied 
from 15" to 57" corresponding to linear scales of 4 - 17 pc. 
At 345 GHz, the beam was slightly broadened in E-W di- 
rection; this is caused by astigmatism, reflecting the lim- 
ited surface accuracy of the SEST antenna (L.-A. Nyman, 
priv. comm.). 

The observations were carried out in a dual beam- 
switching mode (switching frequency 6 Hz) with a beam 



thrc 



of 11'40" in azimuth. The on-source integra- 



In July 1997, a 0.85 mm SIS receiver with T s , 
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tion time of each spectrum varied from 8 minutes for 
12 CO(1-0) to 18.5 hours for C 34 S(3-2). All spectral in- 
tensities were converted to a Tmb scale, correcting for 
main beam efficiencies of 0.74 at 85-100, 0.70 at 100- 
115, 0.67 at 130-150, 0.45 at 220-265GHz, and 0.30 
at 345 GHz (L.B.G. Knee, L.-A. Nyman, A.R. Tieftrunk, 
priv. comm.). Calibration was checked by measurements 
of Orion KL, M17SW, and NGC4945; it was found to 
be accurate to ± 10%. The pointing accuracy, obtained 
from measurements of the SiO masers RDor and UMen, 
was better than 10". 345 GHz data with pointing offsets 
in excess of 5" were ignored in the data analysis. 

3. Results 

Fig. [l] shows maps of the LIRS 36 complex, obtained in 
12 CO J=l-0, 2-1, 3-2, and 13 CO J=l-0 with spacings of 
28", 10", 10", and 28" for the central region, respectively. 
The 1-0 map does not spatially resolve the source. While 
the 2-1 line data are marginally resolving the source in E— 
W direction (full width to half power (FWHP) size « 40", 
corresponding to 12 pc at D — 60 kpc), the emitting region 
is elongated along the N-S axis (FWHP size « 60"; 17 pc 
at D — 60 kpc). The J=3-2 emission shows a similar but 
more compact distribution. Radial velocities are slightly 
above 126 kms -1 toward the central and western parts of 
the complex and 124 - 126 kms -1 at southern, eastern 
and northern offsets. 

A total of 12 molecular species has been observed. We 
have detected nine of these in a total of 18 rotational tran- 
sitions, including isotopic lines. Two tentative detections 
have also been obtained. This demonstrates that it is pos- 
sible to carry out molecular multiline studies, including 
many species, for selected regions of the SMC, extending 
the range of metallicities by more than a factor of two 
below that of the Large Magellanic Cloud (cf. Westerlund 
1990). Spectra and line parameters obtained by gaussian 
fits are displayed in Fig.| and Table 0. The three H 2 CO 
J=2-l K a = 0,1 transitions have been measured, while the 
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Fig. 2. Molecular spectra toward LIRS 36, with the line rest frequency increasing from the upper left to the lower right spectrum. 
Observed nominal peak position: ai 950 = 00 h 44 m 51 s , <5i 95 o = -73°22'30". 



Jk b _,k c — 3o,3-2o,2 transition is only tentatively detected. 
Three rotational transitions of CS, the J=2— 1, 3-2, and 
5-4 lines, have been measured in the main species and two 
in C 34 S. While the C 34 S(3-2) line is apparently detected, 
we find no evidence for the corresponding 2-1 transition. 
Although studies of Galactic star-forming regions (Chin 
et al. 1996) imply that CS (3-2) transitions from rare iso- 
topes are more easily detected than the corresponding 2- 
1 lines, we classify our C 34 S detection as tentative. The 
non-detection of the HCO + and HCN J=3-2 transitions 
is not surprising in view of the relatively high upper lim- 
its obtained in the warm humid weather conditions of the 
Chilean summer. 

Among the rare isotopic species of CO, 13 CO is seen 
but not C 18 0. The tentative detection of C2H (Chin et al. 
1997) is confirmed; the relative line intensities of its two 
observed hyperfine components imply that the emission is 
optically thin (cf. Nyman 1984). CN remains unconfirmed, 
in spite of the higher sensitivity of the data presented here. 

4. Discussion 

LIRS 36 is one of the most prominent far infrared sources 
of the SMC (Schwering & Israel 1989). It is an IRAS point 
source with a cool dust spectrum (Td us t ~ 30 K) and an 
infrared luminosity of 3 x 10 5 Lq (we applied the method 
outlined by Wouterloot & Walmsley 1986; D = 60kpc). 



4-1- The X- factor 

Assuming that the interstellar medium is in virial equi- 
librium at all linear scales and that the line width is rep- 
resenting the cloud's velocity distribution, we can deduce 
the X = Nh 2 /Ico conversion factor from 



M vir = 190 R (Aw 1/2 ) 2 



M„ 



1.36 [AT H2 //co] ™h 2 



i=l 



A, 



and 



M vir =M, 



gas • 



M V i r is the virial mass in solar units (M Q ), M gas is the 
cloud mass derived from Ico (also in solar units), R is the 
cloud radius in pc, Avi/2 is the total linewidth in kms -1 , 
1.36 is the correction to include helium and metals, Ah 2 
is the H2 column density in cm~ 2 , Ico denotes the in- 
tegrated CO (1-0) line intensity in Kkms -1 , toh 2 is t ne 
mass of an H2 molecule, and 7 u .i and Zi ( i (in Kkms -1 ) are 
the integrated intensities of the upper and lower contours 
confining the area A\ (in cm 2 ; n: number of contours). 
The factor 190 in the virial mass equation refers to a 1/r 
density distribution and is consistent, within a factor of 
two, with constant and 1/r 2 density profiles (MacLaren 
et al. 1988). Making use of the beam deconvolved FWHP 
cloud size deduced from the 2-1 spectra (cf. Sect.[|) and 
accounting for all the CO emission from the cloud, we 
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Table 1. Parameters of the Observed Molecular Lines toward LIRS36 



Molecule Frequency Tyg^ r.m.s. b ) i>LSR & v \/2 J ^MB 

& Transition [MHz] [K] [mK] [kms" 1 ] [kms" 1 ] [K km 



C3H2 


J= 2 1, 2-1-0,1 


85338.890 


0.025 


7 


126.3 


3.9 


0.108 ± 0.008 


C 2 H 


N— 1-0 J— 3/2-1/2 F—2 1 


87316.925 


0.040 


8 


126.0 


4.1 


0.180 ± 0.008 




F— 1-0 


87328.624 


0.030 


7 


126.6 


1.8 


0.075 ± 0.008 


HCN 


J— 1-0 


88631.847 


0.026 


5 


126.3 


2.5 


0.106 ± 0.006 


HCO~^ 


J— 1—0 


n 1 cio 

89188.518 


0.101 


17 


126.4 


3 . 1 


0.3by ± 0.018 


HNC 


J— 1-0 


90663.543 


0.016 


5 


126.0 


2.4 


0.038 ± 0.005 


b 


-7=2 1 


96412.982 


< 0.01 


7 






< 0.026 


Ob 


J— 2—1 


nvn n (ico 

y /you.96o 


0.067 


1C) 


126. 2 


2.8 


n noc _i_ n n 1 c 
O.zzo ± O.Ulo 




J— 1—0 


1 nn7o n 1 en 
1U9 f o2.1oU 


< 0.01 


8 






< 0.030 




J — i— u 


1 luzui.ooo 


0.284 


4.5 


126.2 


3.1 


1.U4 It U.UO 


CN 


iv— 1— u J —0/ 2—1/ 2 r—o/2—o/2 


11349U.982 


< 0.01 


12 






< 0.038 


LA J 


J — 1— U 


1 lOi! / 1. ZU4 


2.73 


135 


126. 1 


3.2 


o. kg -I- n 1 A 
y.oy ± U.14 


bU 


J— 43-32 


1001 To C/io 

1381 /W.u4h 


0.051 


1 1 


126.3 


2 . 3 


U. loo it U.U1U 


i±2 OD 




1 /inOQH CIO 

140oo9.51o 


0.081 


lfi 


126. 1 


2.7 


n Oil 1 _i_ r\ m A 
U.z41 it U.U14 


O b 


J — o—Z ' 


1 A Ati~t *7 TAT 

14401 f . 14 f 


0.015 


7 


125.3 


1.9 


O.Olo ± U-UUf 


CH3OH 


J— 3o-2o, A + 


14olUo.zoU 


< 0.01 


10 






< 0.029 


H2 CO 


J= 2 0,2- ln ,l 


I4ob02.yo3 


0.059 


lfi 


126.2 


2.2 


n 1 in 1 n (11 Q 
0.14^ ± U.U13 


cs 


J— 3— 2 


14byb9.04y 


0.065 


14 


126.3 


2.5 


0.179 ± 0.012 


NO 


2n i/2 J=3/2-l/2 F=5/2-3/2 (-+) 


150176.480 


< 0.01 


17 






< 0.048 


H 2 CO 


J= 2 l,l-ll,0 


150498.339 


0.069 


1!) 


126.4 


2.5 


0.167 ± 0.016 


C 3 H 2 


^^ 4 l,4-3o,3 


150851.910 


< 0.02 


20 






< 0.055 


H 2 CO 


^^3o,3~ 2 0,2 d ^ 


218222.191 


0.060 


44 


126.4 


1.7 


0.150 ± 0.034 


13 CO 


J=2-l 


220398.686 


0.565 


48 


126.1 


3.0 


1.85 ± 0.03 


H 2 CO 


^^3i2- 2 i,i 


225697.772 


< 0.04 


74 






< 0.168 


CO 


J=2-l 


230537.990 


3.92 


121 


126.1 


3.2 


13.9 ± 0.2 




J=2-l c ) 


230537.990 


3.02 


112 


126.1 


3.5 


11.7 ± 0.1 


CS 


J=5-4 


244935.606 


0.161 


93 


125.7 


1.6 


0.287 ± 0.067 


HCN 


J=3-2 


265886.432 


< 0.1 


250 






< 0.52 


HCO + 


J=3-2 


267557.625 


< 0.1 


120 






< 0.25 


CO 


J=3-2 


345795.975 


3.23 


570 


125.9 


3.7 


12.7 ± 0.3 




J=3-2 °) 


345795.975 


2.37 


130 


126.1 


3.7 


9.66 ± 0.08 



a) For non-detections, the single channel 3<r noise level is given, divided by the square root of the number of expected line 
channels. 

b) r.m.s noise for an individual channel. 

c) The intensities are integrated from ulsr = 122 to 130 kms" 1 . For the detected lines, errors are one standard deviation. For 
undetected transitions, upper 3<r limits are given. 

d) Tentative detection 

e) The spectrum has been convolved to the beam size of the J— 1-0 transition (~ 43"). 



obtain X LTRS36 = 4.8 x KPcirr 2 (Kkms" 1 )- 1 . This is 
a factor of two larger than the local Galactic disk value 
and leads to a H2 column density (averaged over the map) 
of -/V(H 2 ) ~ 3.7 x 10 21 cm~ 2 , that is consistent with the 
characteristic column density of nearby Galactic molecular 
clouds (e.g. Larson 1981). The cloud follows the correla- 
tions between Avx/ 2 , mass, and size, proposed by Larson 
(1981), and the average number density becomes n(H2) 
« 100 cm- 3 . 

Analysing CO clouds with a wide range of radii, Ru- 
bio et al. (1993a) obtained X SMC « 9 x 10 20 (i?/10pc) 7 
cm -2 (Kkms -1 ) -1 . The correlation between X-factor and 
linear scale R was interpreted in terms of an increased 
rate of photodissociation of CO due to a strong UV radi- 
ation field and a low gas to dust mass ratio and carbon 



abundance. Moreover, Rubio et al. (1993a) suggest a pre- 
dominantly atomic intercloud medium. Our observations 
trace the SMC with a characteristic linear scale of 10 pc 
that is slightly smaller than the size of the molecular com- 
plex; our LIRS 36 X -factor is half of the value obtained by 
Rubio et al. (1993a) from their sample of clouds. 

Specifically for LIRS 36, Rubio et al. (1993b) find a 
virial mass of 5.1 x 1O 4 M0. The mass derived by us 
is 1.6 x 1O 4 M0. The application of a slightly smaller 
lincwidth (3.6 km s" 1 ) than those of Rubio et al. (1993b) 
and Lequeux et al. (1994) (3. Skins" 1 ) can be justified 
by the good agreement of the Avi/ 2 derived from our 
12 CO(1-0), (2-1), (3-2) and 13 CO(1-0) observations (Ta- 
ble [l]) . Since the difference caused by linewidths is negligi- 
ble (« 10%), the discrepancy in M v - lr is caused mainly by 
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the cloud radius i?, 18.6 versus 9.8 pc. The discrepancy in 
R is caused by two effects: Firstly, Rubio et al. (1993a) use 
the CO non-deconvolved 0.4Kkms -1 contour (in units of 
antenna temperature) as the cloud boundary; we use in- 
stead the deconvolved FWHP contour. While Rubio et al. 
(1993a) cloud sizes are based on CO (1-0) data, our cloud 
size is based on the CO (2-1) transition, assuming that 
1-0 and 2-1 emission have a similar extent. Since CO (1— 
0) and (2-1) cloud sizes should be similar for the SMC 
(cf. Lequeux et al. 1994) and since cloud size measure- 
ments based on higher resolution CO (2-1) spectra pro- 
vide more accurate results, our smaller R value should be 
preferred. The determination of R (defined by the FWHP 
contour) from CO (1-0), not (2-1) line spectra, may lead 
to a systematic overestimate of the average A-value for 
clouds with radii R « 10 pc. A small A-factor may im- 
ply that an atomic interclump medium can only play a 
minor role in the molecular cores. Detailed studies of ad- 
ditional SMC cores are needed to demonstrate that small 
scale SMC A-factors are generally as low as is indicated 
by our data. 

4-2. Molecular abundances 
4.2.1. 12 CO versus 13 CO 

Our measured 12 CO and 13 CO line intensities are larger 
than those reported by Rubio et al. (1993b): While the dis- 
crepancy is typically 20%, our 13 CO (2-1) peak line tem- 
perature is higher by a factor of two. This is partially com- 
pensated by the exceptionally large linewidth assigned to 
this line by Rubio et al. (1993b), so that the ratio of in- 
tegrated line intensities only becomes 1.3. Our beam size 
corrected J =2-1/ J =1-0 12 CO and 13 CO line intensity ra- 
tios are 1.21 ± 0.03 and 1.50 ± 0.09 (the la error refers to 
the noise in the individual spectra and does not include 
calibration uncertainties that are given in Sect.§). This 
should be compared with 1.16 and 1.25 from Rubio et al. 
(1993b). Most of the discrepancies in line intensity ratios 
can be explained in terms of the higher sensitivity of our 
A = 1.3 mm data and the smaller 1.3 mm SEST beam ef- 
ficiency adopted by us. 

The integrated line intensity ratios from the presum- 
ably optically thin 13 CO transitions are consistent with a 
density of n(H 2 ) « 10 4 c m - 3 (for T kin « 20 K, applying 
a Large Velocity Gradient (LVG) radiative transfer code; 
see also Lequeux et al. 1994). Since the density is higher 
than the virial density obtained in Sect 



stringent than those obtained from star-forming regions of 
the LMC, the actual value must be larger than the charac- 
teristic Galactic line intensity ratio, J( 13 CO)//(C 18 0) « 
10 (cf. Lada 1976; Johansson et al. 1994). Whether this is 
caused by 'anomalies' in the isotopic abundances relative 
to those of the Galactic disk (cf. Henkel & Mauersberger 
1993), whether it is caused by isotope selective photodis- 
sociation (e. g. van Dishoeck & Black 1988; Fuente et al. 
1993) or whether it is caused by fractionation in a par- 
tially ionized medium (cf. Watson et al. 1976) remains to 
be seen. The small C 18 0/C 17 line intensity ratio in the 
LMC (Johansson et al. 1994) hints at a low 18 abun- 
dance, but we do not know whether this also holds for the 
SMC. 7( 13 CO)//(C 18 0) w 40, as measured toward the 
Galactic Hn region S68 (Bally & Langer 1982), demon- 
trates that fractionation can, in principle, account for the 
observed line intensity ratio anomaly. Assuming that both 
13 CO and C ls O are optically thin and have the same exci- 
tation temperature, that the 13 CO abundance is enhanced 
by the maximum factor permitted by chemical fractiona- 
tion, e 35K / Tk '=, and that the 13 C/ 12 C : ls O/ 16 ratio is w 
7 as in the local Galactic disk (e. g. Wilson & Rood 1994), 
we obtain for 7( 13 CO)//(C 18 0) « 40 with 



a 35 K/T kin 



[7( 13 CO//(C 18 0] / [ 13 C/ 12 C : 18 0/ 16 0] 



a kinetic temperature of « 20 K. The temperature is con- 
sistent with the cloud core model temperatures suggested 
by Lequeux et al. (1994). A large 13 CO abundance en- 
hancement caused by fractionation is plausible in the in- 
terstellar medium of the SMC with its small dust opacities 
and strong UV radiation field. 

4.2.3. HCO+, HCN, and HNC 

The relative J=l-0 intensities of HCO+, HCN, and HNC, 
three molecules with similar rotational constants and elec- 
tric dipole moments, also follow the trend obtained toward 
star-forming regions of the LMC. For the J =1-0 transi- 
tion, 7(HCO+)//(HCN) and 7(HCN)/J(HNC) > 1. This 
can be interpreted in terms of warm molecular gas, cou- 
pled with an intense ionization flux from supernovae, and 
with HCO + arising from a larger volume than HCN (cf. 
Johansson et al. 1994; Chin et al. 1997). 



4.1 



there must 4.2.4. CN 



be small scale structure that is not resolved by our <; 15" 
beam. 

4.2.2. 13 CO versus C 18 

The non-detection of C 18 is consistent with the large 
13 CO/C 18 line intensity ratios observed toward the 
LMC (Johansson et al. 1994; Chin et al. 1997). We find 
/( 13 CO)//(C 18 0) > 35 (a 3a limit). While the limit is less 



The most unexpected result of our line survey is the ab- 
sence of a detectable CN signal. CS J=2-l and 3-2 lines 
tend to be weaker than those of CN in most Galactic and 
extragalactic sources (cf. Henkel et al. 1988, 1990, 1993; 
Mauersberger & Henkel 1989; Nyman & Millar 1989; Zi- 
urys et al. 1989; Nyman et al. 1993). Toward LIRS 36, CN 
remains undetected but we may have seen a rare isotopic 
species of CS. 
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4.2.5. CS and the cloud density 

At least three of the five observed CS transitions were de- 
tected (see Table [l] and Fig. ||). In the case that the CS 
emitting region is as extended as the CO emitting region, 
we have to apply the CO 1.3 mm line intensity correction 
factor, 0.7, for CS J— 5-4. This factor can be extracted 
from Table Q and was also recommended by Rubio et al. 
(1993b) to account for the small 1.3 mm beam size when 
comparing 1.3 with 3 mm data. Interpolating, 0.85 is then 
the appropriate correction factor for our A = 2 mm data. 
If CS arises instead from a point source, the beam size 
corrections become 0.25 at 1.3 mm and 0.5 at 2 mm wave- 
length. Since CS requires higher densities than CO to be- 
come detectable and since CS optical depths tend to be 
smaller (cf. Bohlin et al. 1978; Linke & Goldsmith 1980; 
Larson 1981; Bachiller & Cernicharo 1986), actual correc- 
tion factors will be 0.25 < /1.3 < 0.7 at 1.3 mm and 0.5 
< /2 < 0.85 at 2 mm. Applying an LVG model describ- 
ing a spherical or planparallel cloud with uniform density 
and temperature (for details, see Mauersberger & Henkel 
1989), we then find that in the optically thin case with Tkm 
= 20 K (cf. Lequeux et al. 1994) the number density be- 
comes 10 6 cm~ 3 < n(H2) < 10 66 cm~ 3 . Also accounting 
for the standard deviation of the integrated CS (5-4) line 
intensity (Table|l|), we find n (H 2 ) = 10 5 - 9 " 6 - 9 cm" 3 . With 
a high 1.3 mm beam efficiency of 0.6 (0.45 is used through- 
out this paper; see Sect.0) we would get instead n (H2) = 



10 



5.7-6.4 



cm . In order to fully explore the range of rele- 



vant kinetic temperatures, we have also calculated densi- 
ties for Tkin = 100 K, a very high value (see e. g. Mauers- 
berger et al. 1990; Lequeux et al. 1994). Accounting again 
for the standard deviation in the integrated CS (5-4) line 
intensity, we now find n(H2) = 10 5 2-5 7 cm~ 3 . With all 
uncertainties considered, we thus conclude that densities 
must be high, of order 10 5 to 10 7 cm~ 3 . This is consis- 
tent with the density derived from H2CO if the Jko,,Kc 
= 3o,3 — 2o,2 transition is detected and if th 2 co,j=2-i 
1. The densities are much higher than those derived from 
12 CO in Sect. Ill] and from 13 CO in Sect.|4.2.l| and are 



mainly based on a comparison of 2-1 with 5-4 CS lines. 
Because of their relatively large separation in excitation 
energy, they are much more sensitive density tracers than 
e. g. the 2-1 and 3-2 transitions. 

So far, we have assumed that CS lines are optically 
thin, which was suggested by Johansson et al. (1994) for 
N 159, an LMC H 11 region that exhibits stronger line emis- 
sion than LIRS36. With a C 32 S/C 34 S abundance ratio of 
w 20 in the disk of the Milky Way however (see Chin et al. 
1996), a comparison of C 32 S and C 34 S (3-2) line intensi- 
ties (Table[|) leads to an optical depth of ~ 5 for the C S 
(3-2) line. There are hence three possibilities: (1) The CS 
lines are moderately optically thick (r ^ 5). With n(H2) 
> 10 4 ' 5 cm~ 3 , this does not imply a significant reduction 
of the densities estimated above. (2) The CS lines are op- 
tically thin and the 32 S/ 34 S ratio drastically deviates from 



that of the Galactic disk. This leads to a 32 S/ 34 S isotope 
ratio of ~ 5, implying that in the SMC massive stars and 
type la supernovae provide an isotopic mixture that is 
quite different from that of the more metal rich disk of 
the Milky Way. (3) Our tentative C 34 S (3-2) detection, 
although appearing to be quite convincing, is not real. 
Another independent measurement of C 34 S is required to 
discriminate between these possibilities. 



4.2.6. Column densities 



Having estimated the molecular densities in Sects. 4.1 
4.2.1, and 4.2.5 and having a quantitative idea of the ki- 
netic temperature distribution (Lequeux et al. 1994), we 
can now also determine the column densities for various 
molecular species. Both temperature and density vary in- 
side the cloud and each molecule and transition is trac- 
ing gas at different temperature and density; error limits 
are thus difficult to assess. Nevertheless, our data allow 
a rough estimate of molecular abundances in the opti- 
cally thin limit. This is the appropriate approach, firstly 
because we do not know the optical depth of most lines 
and secondly because Johansson et al. (1994) find HCO + , 
HCN, and CS to be (almost) optically thin in their mm- 
wave transitions toward the prominent LMC star-forming 
region N 159. Assumed kinetic temperatures, densities, 
and resulting absolute and fractional (relative to H2) col- 
umn densities for a ss 50" (15 pc) beam size are given in 
Table I 

A comparison with Galactic disk values directly shows 
that the LIRS36 fractional abundances, [X]/[H2], are 
much smaller (e. g. Blake et al. 1987): For CO, CS, HCO+, 
and H2CO, underabundances relative to the Orion ex- 
tended ridge or to the prototypical dark cloud TMC-1 are 
w 10. For HCN and HNC and probably for CN as well, 
underabundances reach two orders of magnitude and are 
thus even surpassing those reported by Johansson (1997) 
for the LMC star-forming region N 159. It is quite re- 
markable that CN is not enhanced relative to HCN and 
HNC as observed in Galactic disk photodissociation re- 
gions with high UV fluxes (Fuente et al. 1993; Greaves 
& Church 1996). We may suspect that the low LIRS36 
fractional abundances of HCN, HNC, and CN are related 
to a small nitrogen abundance. With average values of 
[N/H] = — 1.1 and [C/H] = —1.3 in star- forming regions of 
the SMC (Westerlund 1990), however, nitrogen is not the 
most depleted among the CNO elements. It thus remains 
to be seen whether the HCN, HNC, and CN fractional 
abundances are related to an exceptionally small nitrogen 
abundance in LIRS 36 or whether they are caused by a 
chemical process that is found in other molecular SMC 
cores as well. 
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Fig. 3. Chemical abundances of selected molecules: SMC molecular abundances for densities n (H2) of (a) 10 4 , (b) 10 5 , and (c) 
10 6 cm -3 as a function of chemical timecscale, using elemental abundances from the 'S2' model of Millar & Herbst (1990); (d) 
'Galactic' model (model 'G' of Millar & Herbst 1990). 



4.2.7. Chemical model calculations 

We have used the models described by Chin et al. (1997) 
to investigate the chemical evolution of LIRS 36 using el- 
emental abundances of C, N, and O from the S2 model of 
Millar & Herbst (1990). Because observed emission clearly 
arises in regions having different densities and because 
we have no information on the density profile within the 
molecular cloud nor on the the size of the high density 
CS clump, we have considered a simple constant density 
model with densities in the range n (H2) = 10 4 — 10 6 cm~ 3 
and Tkin = 20 K. For comparison, we have also calculated 
a 'Galactic' model (model 'G' of Millar & Herbst 1990) for 
^(H^) = 10 4 cm~ 3 (Fig. ||). In addition to a lower metal- 
licity, the SMC models have a smaller dust-to-gas ratio, 
and hence H2 formation rate, by a ratio of 6.5 and a UV 



intensity larger by a factor of 4, compared to Galactic 
values. There are conflicting results on the appropriate 
cosmic-ray ionisation rate. Chi & Wolfendale (1993) have 
argued that, based on gamma-ray observations, the rate 
in the SMC is at most 11% of the Galactic value. How- 
ever, observations of HCO+ in the LMC and SMC (Chin 
et al. 1997) indicate that the rate is larger than that im- 
plied by Chi & Wolfendale. In these models we use a rate 
of 5.2 x 10 _17 s _1 , three times larger than the Galactic 
value. The main impact of this rate is to increase the ion- 
isation fraction and to reduce the chemical time-scales. 
Finally, we have used a visual extinction of 10 mag. in 
calculating photorates. While this is obviously too large 
for the molecular envelope of LIRS 36, the central region, 
in which densities greater than 10 5 cm -3 are derived from 
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Table 2. Molecular column densities 



Molecule 


T kin 

[K] 


n(H 2 ) 
[cm" 3 ] 


Column density 
[cm- 2 ] 


Fractional abundance 


H 2 a ) 






3 X 


10 21 




1 


12 co b) 






3 x 


10 16 


1 


x 10- 5 


13 GO b) 


20 


1Q 3.9 


7 x 


10 14 


2 


x to -7 


CN c) 






< 2 x 


10 12 


< 7 


x 10- 10 


CS d > 


20 


10 65 


1 X 


10 12 


3 


x 10- 10 


HCO+ c ~> 


30 


10 5.0 


7 X 


10 11 


2 


x 10- 10 


HCN f > 


20 


10 6 5 


2 x 


10 11 


7 


x 10- 11 


HNC g) 


20 


10 5.5 


9 x 


10 10 


3 


x 10- 11 


H 2 CO h) 


20 


10 5.5 


9 x 


10 11 


3 


X to" 10 



Colu mn density 
Sect. 



b) 



4.1) 



estimated from the virial mass (see 



d) 
e) 



f) 



The column density has been derived from an LVG simula- 
tion, reproducing the beam size corrected 13 CO J =1-0, 2- 
1, and 3-2 line intensities (Tableland Sect. 1.2.1 ) and ap- 
plying a velocity gradient of Awi/2/2-R = 0.27 kms -1 pc _1 . 
With a 12 CO/ 13 CO abundance ratio of 50 (Johansson et al. 
1994) the CO column density becomes 3.5 x 10 16 cm" 2 . The 
fractional abundance is [CO]/[H 2 ] = 10" 5 . This is « 20% 
of the value in Galactic clouds (e. g. Wilson et al. 1986), re- 
flecting the low carbon abundance in star-forming regions 
of the SMC ([C/H] » -1.3; Westerlund 1990). 
CN excitation temperatures are small (< 10 K) at low den- 
sities (n(H 2 ) < 10 4 cm -3 ), while abundances decrease at 
higher densities (e. g. Turner & Gammon 1975; Churchwell 
1980; Churchwell & Bieging 1982, 1983; Crutcher et al. 
1984; Henkel et al. 1988; Johansson et al. 1994). With an 
excitation temperature of 10 K and accounting for the fact 
that the LTE intensity of the N=l-0 J=3/2-l/2 F=5/2- 
3/2 CN line only comprises one third of the total N=l-0 
line intensity under optically thin conditions (Skatrud et 
al. 1983), we find the given beam averag ed 3a upper CN 
column density limit. See also Sect. 
See Sect.|4.2.5. 



1.2 



h) 



The column density is derived in the same way as for HCN 
(see footnote f). Different densities and temperatures were 
chosen to account for HCO + depletion at highest densities 
and for the suspected large spatial extent of the HCO + 
emission region (e.g. Chin et al. 1997). 
HCN is suspected to emit mainly from the coolest and dens- 
est parts of the cloud (cf. Johansson et al. 1994; Chin et 
al. 1997), so that T kin = 20 K and n(H 2 ) = 10 6 5 cm -3 
are appropriate. The column density was derived with an 
LVG code (see also footnote b). The HCN (3-2) upper limit 
(Tableg) is not small enough to provide a stringent density 
limit. 

The column density is derived in the same way as for HCN 
(see footnote f). Since [HCN] /[HNC] abundance ratios tend 
to drop with rising density (Schilke et al. 1992), n (H 2 ) = 
10 5 - 5 cm -3 is justified. 

The column density was derived from an LVG simulation 
with J=3-2 line intensities being sli ghtly weaker (by ~ 
30%) than the 2-1 lines (see also Sect. 4.2.5[). 



CS, could easily have an extinction large enough to ex- 
clude photoeffects even with a low dust-to-gas ratio. 

The results of our four model calculations are sum- 
marised in Fig. U and indicate that the abundance of CN 
is sensitive to the O2 abundance since this species destroys 
CN with a fast, measured rate coefficient of 2.4 x 10~ n - 
(r kin /300K)-°- 6 cm 3 s- 1 down to 20 K (Sims et al. 1994). 
The fractional abundance of O2 increases as one goes from 
the Galactic to the SMC model at n (H 2 ) = 10 4 cm -3 and 
as density increases because of the decreasing abundance 
of C atoms which destroy O2. In model S2 there is more 
free oxygen remaining after all available carbon has been 
tied up in CO so that the importance of the destruction 
of CN by O2 increases. Note that the O/C ratio, where X 
represents the elemental abundance of species X, varies 
between 2.4 (model G) and 10.5 (model S2) while the dif- 
ference between the elemental abundances of O and C is 
2.1 x 10_~ 4 (model G) and 5.7 x 10 -5 (model S2), and the 
ratio (O — C)/C, increases from 1.4 (model G) to 10.4 
(model S2). These values show that in model S2 there is a 
lot of oxygen remaining after all the carbon has been pro- 
cessed into CO and thus the importance of CN destruction 
by 2 is more important in the SMC models. 

The fractional abundance of CN decreases by about 
a factor of 3 as one goes from model G to S2 and by 
a further factor of 100 as the density increases from 10 4 
to 10 6 cm~ 3 . As a result, the abundance of CN becomes 
less than or comparable with those of HCN and HNC at 
high density. On the other hand, CS is not destroyed by 
O2 or O atoms and its abundance changes by less than a 
factor of two over all the models investigated. Note that 
the maximum fractional abundance of CN at 10 6 cm -3 is 
only 4 x 10~ 10 for chemical timescales < 10 5 yrs; for longer 
times, the abundance falls rapidly to its steady-state value 
of 4 x 10 -14 . 

5. Conclusions 

Having made a mm- wave molecular survey of the LIRS 36 
star- forming region in the SMC, we obtain the following 
main results: 

(1) We have detected CO, CS, SO, C 2 H, HCO+, HCN, 
HNC, H 2 CO, and C 3 H 2 . NO, CH 3 OH, and, surpris- 
ingly, CN were not detected. Among the rare isotopic 
species, we have seen 13 CO and tentatively C 34 S, but 
not C 18 0. 

(2) For a characteristical scale length of » 10 pc and as- 
suming virial equilibrium, the X = N{S.2)/Ico con- 
version factor is with w 4.8 x 10 20 cm -2 (Kkms -1 ) -1 
a factor of two larger than the local Galactic disk X- 
factor. SMC JT-factors given by Rubio et al. (1993a) 
for this scale length may have to be reduced by half an 
order of magnitude. 

(3) Density estimates range from n(H2) ~ 100 cm -3 (de- 
duced from the virial masss and the spatial extent of 
the CO emission) over 10 4 cm -3 (from 13 CO) to 10 5 - 
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10 7 cm -3 (from CS). The observed CS transitions pro- 
vide strong evidence for the presence of a very dense 
core with a density likely surpassing 10 6 cm~ 3 . 

(4) As in the LMC cloud cores, the 7( 13 CO)//(C 18 0) line 
intensity ratio is larger than the usual values encoun- 
tered in the Galactic disk. Whether this is caused 
by isotopic abundance anomalies, by isotope selec- 
tive photodissociation, or by chemical fractionation re- 
mains an open question. 

(5) 7(HCO+)/J(HCN) and J(HCN)/7(HNC) line inten- 
sity ratios are > 1, consistent with molecular emission 
from a warm (Tki n > 10 K) molecular environment ex- 
posed to a high ionizing flux. 

(6) Fractional (relative) abundances of CO, CS, HCO + , 
and H2CO are an order of magnitude below those 
of the Galactic disk. HCN, HNC, and likely also CN 
are even underabundant by two orders of magnitude. 
Whether this is reflecting a particularly small nitrogen 
abundance in LIRS 36 or whether this is a common 
chemical peculiarity of other SMC cores as well re- 
mains open. The non-detection of CN can be explained 
in terms of the high cloud density and a (relatively) 
high fractional abundance of O2 that is destroying CN 
but not CS. 
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